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Abstract: We calculate the drag force experienced by an infinitely massive quark prop- 
agating at constant velocity through an anisotropic, strongly coupled = 4 plasma by 
means of its gravity dual. We find that the gluon cloud trailing behind the quark is gen- 
erally misaligned with the quark velocity, and that the latter is also misaligned with the 
force. The drag coefficient ^ can be larger or smaller than the corresponding isotropic 
value depending on the velocity and the direction of motion. In the ultra-relativistic limit 
we find that generically fi (x p. We discuss the conditions under which this behaviour may 
extend to more general situations. 
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1 Introduction 

A remarkable conclusion from the experiments at the Relativistic Heavy Ion Collider 
(RHIC) [1, 2] and at the Large Hadron Collider (LHC) [3] is that the quark-gluon plasma 
(QGP) does not behave as a weakly coupled gas of quarks and gluons, but rather as a 
strongly coupled fluid [4, 5]. This renders perturbative methods inapplicable in general. 
The lattice formulation of Quantum Chromodynamics (QCD) is also of limited utility, since 
for example it is not well suited for studying real-time phenomena. This has provided a 
strong motivation for understanding the dynamics of strongly coupled non-Abelian plas- 
mas through the gauge/string duality [6-8] (see [9] for a recent review of applications to 



For a period of time immediately after the collision, the system thus created is 
far from equilibrium. After a time Ti^^ > "^out the system becomes locally isotropic and 
a standard hydrodynamic description becomes applicable. It has been proposed than an 
intrinsically anisotropic hydrodynamical description can be used to describe the system at 
intermediate times r^nt < r < Ti^o [10-18]. In this phase the plasma is assumed to have 
significantly unequal pressures in the longitudinal and transverse directions. The standard 
hydrodynamic description is a derivative expansion around equal pressures, and therefore 
it is not applicable in this regime. In contrast, the intrinsically anisotropic hydrodynamical 
description is a derivative expansion around an anisotropic state, and hence in this case the 
requirement that derivative corrections be small does not imply small pressure differences. 
In a real collision the degree of anisotropy will decrease with time, but for some purposes 
it is a good approximation to take it to be constant over an appropriate time scale. 



the QGP). 
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Motivated by these considerations, in this paper we will investigate the effect of an 
intrinsic anisotropy on the drag force felt by a quark propagating through a strongly coupled 
plasma. For this purpose we will examine a string moving in a gravity solution [19, 20] dual 
to an anisotropic = 4 super Yang-Mills plasma. As we will review below, the plasma 
is held in anisotropic equilibrium by an external force. The gravity solution possesses an 
anisotropic horizon, it is completely regular on and outside the horizon, and it is solidly 
embedded in type IIB string theory. For these reasons it provides an ideal toy model in 
which questions about anisotropic effects at strong coupling can be addressed from first 
principles. 

Following the original calculations [28, 29] of the drag coefficient, the closely related 
diffusion coefficient was obtained independently in [21]. These seminal papers have been 
generalized and elaborated on in a vast number of subsequent contributions [22] , including 
in particular comparisons with the corresponding weakly-coupled results [23], as well as 
extensive analyses of the energy-momentum tensor which provide a detailed picture of the 
directionality of energy flow away from the moving quark [24]. Examples of holographic 
studies of the drag force in the presence of anisotropies and/or inhomogeneities include 
[25, 26]. 



2 The gravity solution 

The type IIB supergravity solution of [19, 20] in the string frame takes the form 

ds^ = ^ (-TB + dx^ + dy' + ndz' + ^)+ L^e'^Unl (2.1) 
X = az, ^ = Hu), (2.2) 

where x 4' ^-^s the axion and the dilaton, respectively, and (t, x, y, z) are the gauge 
theory coordinates. Since there is rotational invariance in the xy-directions, we will refer 
to these as the transverse directions, and to z as the longitudinal direction. J^, B and % 
are functions of the holographic radial coordinate u that were determined numerically in 
[19, 20]. Their form for two values of a/T is plotted in Fig. 1. The horizon lies at n = Uh, 
where = 0, and the boundary at n = 0, where J-=B = 'H = (f) = l. The metric near 
the boundary asymptotes to AdS^ x S^. Note that the axion is linear in the z-coordinate. 
The proportionality constant a has dimensions of mass and is a measure of the anisotropy. 
The axion profile is dual in the gauge theory to a position-dependent theta parameter of 
the form 9 (x z. This acts as an isotropy-breaking external source that forces the system 
into an anisotropic equilibrium state. 

If a = then the solution reduces to the isotropic black D3-brane solution dual to the 
isotropic A/" = 4 theory at finite temperature. In this case 

1 

B = n = l, X = (^ = 0, T=l-^, u^ = — (2.3) 
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Figure 1. Metric functions for a/T ~ 4.4 (left) and a/T ~ 86 (right). 




Figure 2. Log-fog pfot of the entropy density as a function of a/T, with Si^o defined as in eqn. (2.4). 
The dashed blue line is a straight line with slope 1/3. 

and the entropy density takes the form 

s,.„ = y7V2r3. (2.4) 

Fig. 2 shows the entropy density of the anisotropic plasma as a function of the dimen- 
sionless ratio a/T, normalized to the entropy density of the isotropic plasma at the same 
temperature. At small a/T the entropy density scales as in the isotropic case, whereas at 
large a/T it scales as [19, 20, 27] 

. = c.„,iV2„l/3r8/3 , (2.5) 

where c^nt is a constant that can be determined numerically. 

A feature of the solution (2.2) that played an important role in the analysis of [19, 20] 
is the presence of a conformal anomaly. Its origin lies in the fact that diffeomorphism 
invariance in the radial direction u gets broken in the process of renormalization of the on- 
shell supergravity action. In the gauge theory this means that scale invariance is broken 
by the renormalization process. One manifestation of the anomaly is the fact that, unlike 
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the entropy density, other thermodynamic quantities do not depend solely on the ratio a/T 
but on a and T separately. Fortunately, this will not be the case for our drag force, which 
will take the form -F(a, T) = T^/ with / a function of the ratio a/T alone. The reason for 
this is that no regularization procedure is necessary for the computation of the drag force, 
and thus diffeomorphism invariance is preserved. We will also verify this analytically in 
certain limits, and numerically for general values of a and T. 



3 Drag force 

Extending the isotropic analysis of Refs. [28, 29], in this section we will consider the drag 
force acting on a quark moving at constant velocity through the anisotropic N = A plasma 
described by (2.2). A simple model for this system is described by the equation of motion 

f=-MP+F, (3.1) 

where p is the quark's momentum, /i is a drag coefficient, and F is an external force. The 
necessary force to keep a steady motion \s F = ^p. An observation that will be important 
for us is that, in the case of an anisotropic medium, the drag coefficient is not just a number 
but a matrix. In our case we will see that this matrix is diagonal, fi = diag{fix, ^J'y, fJ-z) with 
= Aij/ 7^ /U^. Thus we should expect that the force and the momentum or the velocity 
of the quark will not be aligned in general, and indeed our calculations will reproduce 
this feature. We will also see that, unlike in [28, 29], the drag coefficient is momentum- 
dependent. 

On the gravity side the quark is described by a string propagating in the background 
(2.2). The string action is 

where g is the induced worldsheet metric. With the factor from the spacetime metric 
the Lagrangian scales as /2-Ka' = \/A/27r. We will set this factor to one in intermedi- 
ate expressions, and we will reinstate it at the end. Denoting the spacetime coordinates 
collectively by X*^, the flow of spacetime momentum 11 a/ along the string is given by 

Physically, one can imagine that the external force on the quark needed to sustain steady 
motion may be exerted by attaching the endpoint of the string to a D7-brane and turning 
on a constant electric field Fmn = d^M^N] on the brane. In other words, we add to the 
action (3.2) the boundary term 

5bd.. = -/ drANdrX^ = -l [ drFMNX^drX"" . (3.4) 

^ JdT. 
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Demanding that the boundary term arising from variation of the total action S + S'bdry 
vanish yields the boundary condition 

UM + FMNdrX^\g^=0. (3.5) 



We will now specify to the case of a quark moving steadily through the plasma. The 
string will not move along the sphere directions, so this part of the metric will play no role 
in the following. Also, given the rotational symmetry in the xy-directions, we will assume 
that y = 0. We fix the static gauge by identifying {t,u) = {a^,a^) and consider a string 
embedding of the form 



x{t,u) — (^vt + x{u)^ sinip , 
z{t, n) — )• ivt + z{u) ) cos (f , 



(3.6) 
(3.7) 

corresponding to a quark moving with velocity v in the xz-plane at an angle ip with the 
z-axis. Under these circumstances the Lagrangian takes the form 



2^/2 



+ sin^ If {BT^x 



% COS^ if 



Tv^ix' — z'Y sin^ if) 



1/2 



(3.8) 

The rates at which energy and momentum flow down the string towards the horizon are 
then 



1 



1 

£^ 
1 



BTv X sin^ Lp + Hz cos^ 99 
BT X + Hv^iz — x') cos^ f 



n 



BT z + V {x — z) sin if 



simp , 



cos if . 



(3.9) 



where ' denotes differentiation with respect to ti, and the boundary conditions (3.5) become 



-lit = FxV sin if + FzV cos 93 . 



(3.10) 



where {Fx,Fz) denote the components of the external force (the electric field). The first two 
equations are the statement that the external force exactly compensates for the momentum 
lost by the quark into the medium. The third equation is identically satisfied by virtue of 
(3.9), and it expresses the fact that the work done by the external force precisely equals 
the rate at which the quark deposits energy into the medium. As we will see below, the 
energy and the momentum flow from the boundary to the horizon (i.e. na;,!!^ and —Tit) 
are positive provided the string trails behind the quark (i.e. if x\z' are negative), as we 
would expect on physical grounds. This can be easily seen by inspection in the simple cases 
of motion along the z-direction (99 = 0), for which 



BFUz' 



u^JB-^ + BFnz'^ 



-Ut = Ii,v. 



0. 



(3.11) 
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and of motion along the X-direction {ip = '7r/2), for which 

BFx' 



0. 



(3.12) 



We will now determine the string profile and the corresponding values of the energy 
and momentum flows for arbitrary ip. The first observation is that, generically, the string 
does not trail 'below' its endpoint's trajectory. In other words, x{u) ^ z{u). Indeed, if 
x{u) = z(u) then the ratio of the momenta would be given by 



tanip 



(3.13) 



which would be a contradiction because the left-hand side is constant whereas the right- 
hand side is not. In order to determine the correct string profile we invert the relations 
(3.9) to find 



where 

Nx = — na;(SJ^csc — sin 99) + n^u^ cos (/? , 

Nz = —Tiz{BF sec if — T-Lv^ cos if) + YixT~Lv'^ sin ip 
BT esc if sec (f 



'j^VMVNxNz-D 



D 



Hx^zU^ — T~iv^ COS (f sin ip 



BT — {H cos^ 93 + sin^ y?) 



(3.14) 

(3.15) 
(3.16) 

(3.17) 



The factor NxNz — D inside the square root in the denominator of (3.14) is positive at 
the boundary, where B,J-',T-L^1 and u — t- 0, and also at the horizon, where — t- 0, and 
generically it becomes negative in some region in between. In other words, it vanishes at 
two different values of u between the boundary and the horizon. To see this, consider the 
last factor in square brackets in (3.17). BJ^ ((H) is monotonically decreasing (increasing) 
from the boundary to the horizon, so this factor is positive at the boundary and negative 
at the horizon. Therefore there exists a critical value Uc in between such that 



BcJ^c — V (Tic cos if + sin 93) = , 
where Be = B{uc), etc. At this point L> = and 

NxNz\^^ = -v'^ {Uc^x cos if - Hz sm iff , 

which is negative unless the momenta are related through 

tan 93 



-He 



(3.18) 



(3.19) 



(3.20) 



If this condition is not satisfied then NxNz — D is negative in some interval ui < Uc < U2 
and vanishes at u = ui and at u = U2- This type of solutions correspond to strings with 
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two endpoints at the boundary. Here we wish to study isolated quarks, which are described 
by strings that extend ah the way from the boundary to the horizon, so we must require 
that NxNz — D is non-negative for all < u < Uh- This is satisfied if and only if (3.20) 
holds and if the two zeros of D coincide with one another, i.e. if the first square bracket 
in (3.17) also vanishes at u = Uc- The latter condition, together with (3.20), allows us to 
solve for the two momenta independently with the final result: 

„ V sin 03 ^ „ , V cos oj 

= — 2^ , iiz = nc — 2^ . 3.21 

Under these circumstances the denominator in (3.14) is always real and positive except at 
Uc, where it vanishes. At this point the numerators also vanish and the functions 
smooth and negative for all < n < Uh provided in (3.14) we choose the positive sign for 
u < Uc and the negative sign for u > Uc- 

In summary, we have obtained the force F = {Ilx,Ilz) that must be exerted on the 
quark in order to maintain its stationary motion, 



^/A V 
2tt u 



2 (sin cos (/?) , (3.22) 



in terms of the quark's velocity v = f (sin cos (/?). (In this equation we have reinstated 
the factor L^/27rQ'.) The external force F is equal to minus the drag force exerted on the 
quark by the plasma, but in a slight abuse of language we will refer to F itself as the drag 
force. Note that v and F are not aligned with one another except in the isotropic case, for 
which %c = 1, or in the cases in which the velocity is aligned with one of the axes, in which 
(f = 0,7r/2. Note also the force depends on the velocity both through the explicit factors 
of V and (p in eqn. (3.22) and implicitly through the value of T-Lc, which is a solution of the 
iJ-dependent equation (3.18). 

Substituting the result (3.21) in (3.14) we obtain the form of the string profile as a 
function of the velocity. The projection of this profile on the xz-plane has tangent vector 
= {tx^Tz) = {x' sin ip, z' cos if). The angle ipr between this vector and the z-axis is 

tanv?^ = — = etanv9, e = 1 + — — „ ■ (3.23) 

Tz dJ- He - OcJ-crL 

At the horizon we have = and thus e = 1, which means that deep in the infrared the 
string aligns itself with the velocity. However, near the boundary — )• 1 and thus 

This is different from unity for generic (p, and so the string does not align itself with the 
velocity except if ip = or ip = tt/2. In these two special cases the entire string profile (not 
just the infrared part) aligns itself with the corresponding z- or the x-axis, respectively, 
because e remains finite in these limits whereas tanp — t- 0, cxd, respectively. 
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The formulas above reduce to the correct expressions in the isotropic hmit (2.3). In 
this case eqn. (3.18) yields 



ul = uUl-v^ = ^^l^ff- (3-25) 
and the force (3.22) becomes 

FUT) = FUT){smip,cosip) (3.26) 

with 

FUT) = lVXT'-jJL=, (3.27) 
as in [28, 29]. For later purposes it is useful to rewrite this result as 

a/a 5^/3 y 

FUs) = ^ , (3.28) 

^ ^ (27r)i/3iv53 ^ ' 

in terms of the entropy density (2.4) of the isotropic = 4 plasma. 



2 



4 Results 



With the groundwork above in place, we can now proceed to state our results. Since for 
general a the metric functions in (2.2) are only known numerically, we have numerically 
determined the drag force as a function of the magnitude of the quark velocity v, of its 
direction (/?, and of the anisotropy a measured in units of the temperature T or in units of 
the entropy density s. The reason for working with both a/T and a/ s^/^ is that we wish to 
compare the drag force in the anisotropic plasma to that in the isotropic plasma, and this 
can be done at least in two different ways: the two plasmas can be taken to have the same 
temperatures but different entropy densities, or the same entropy densities but different 
temperatures. 

The drag force F{v^ a/T) in units of the isotropic drag force in a plasma at the same 
temperature is shown in Fig. 3. The drag force F{v, ^p, a/ s^^^) in units of the isotropic drag 
force in a plasma at the same entropy density is shown in Fig. 4. With a few exceptions, 
the results are qualitatively similar. In both cases we see that the anisotropic drag is larger 
than the isotropic drag except in a region near the x-axis. This region is more clearly shown 
in Fig. 5: the curves in that figure are the intersections between the two surfaces shown in 
each of the corresponding 3D plot in Figs. 3 or Figs. 4. Considering that the value of a/T 
varies by a factor of 62 between the top and the bottom curves in Fig. 5, we see that the 
region in question depends relatively mildly on the magnitude of the anisotropy. 

For motion along the longitudinal z-direction, the anisotropic drag is greater than the 
isotropic drag for any value of v. For any direction of motion / 7r/2, the ratio FanisoZ-^iso 
diverges as l/Vl — v'^ in the ultra-relativistic limit ?; — )■ 1 irrespectively of whether the 
comparison is made at the same temperature or at the same entropy density, as we prove 
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Figure 3. Drag force as a function of the quark velocity {v^^v^) = u(sin cos (p) for a quark 
moving through an anisotropic plasma with a/T = 1.38(a), 4.41(b), 12.2(c), 86(d). F is plotted in 
the appropriate units to facilitate comparison with the isotropic result (3.27) for a plasma at the 
same temperature. 



analytically in Appendix A. In other words, for motion not perfectly aligned with the 
transverse x-direction, the anisotropic drag becomes arbitrarily larger than the isotropic 
one as the ultra-relativistic limit is approached. This is most clearly illustrated in Fig. 6, 
which shows constant-(/9 slices of the (c) and (d) plots in Figs. 3 and Figs. 4. We will come 
back to this result in Sec. 5. 

The ratio -Famso/-^iso is always finite for motion along the transverse x-direction. (In- 
cidentally, this implies that the limits (/?—)■ vr/ 2 and v — )• 1 do not commute.) In this case 
we must distinguish between the comparisons at equal temperature or at equal entropy 
density. In the first case, our numerical results indicate that the anisotropic drag is smaller 
than the isotropic one for < v < Vc and larger than the isotropic one for Vc < v < 1, 
and we have confirmed this analytically in the limits of small and large anisotropics (see 
the Appendices). The velocity Vc at which the transition takes place is Vc — 0.9 for small 
anisotropics and it approaches 1 as the anisotropy increases. 
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Figure 4. Drag force as a function of the quark velocity {v^^v^) = u(sin cos (p) for a quark 
moving through an anisotropic plasma with a/T = 1.38(a), 4.41(b), 12.2(c), 86(d). F is plotted in 
the appropriate units to facilitate comparison with the isotropic result (3.28) for a plasma at the 
same entropy density. 



In the second case our numerical results indicate that the anisotropic drag is smaller 
than the isotropic one for all v G [0, 1] provided a/s^^^ is small enough. In the opposite limit, 
a/s^^^ S> 1, the anisotropic drag stays smaller than the isotropic one for small velocities 
and becomes larger above some critical velocity. We have confirmed this analytically in 
the Appendices. For a fixed v, the angle with respect to the z-direction beyond which the 
anisotropic drag may become smaller than the isotropic drag is shown in the constant-f 
slices of Fig. 7. 

The dependence of the drag force on the anisotropy for fixed velocity is most clearly 
seen in Figs. 8 and Figs. 9, where the ratio F/Fi^o is plotted for several values of v and ip. 

In order to illustrate the geometric properties of the string solution, in Fig. 10 we have 
plotted the projection of the string profile onto the gauge theory directions. As anticipated, 
we see that the string curves in the xz-plane and (unless 93 = or 7r/2) only aligns itself 
with the velocity in the far infrared, i.e. at large u. The misalignment between the velocity 
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Figure 5. Values of the velocity at which the drag in an anisotropic plasma with (from top to 
bottom) a/T = 1.38,4.41, 12.2,86 equals the drag in an isotropic plasma at the same temperature 
(left) or at the same entropy density (right). For a given value of a/T, the anisotropic drag is larger 
(smaller) than the isotropic drag above (below) the corresponding curve. 

V, the drag force F, and the tangent to the string profile at the string's endpoint r are 
shown in Fig. 11. We see that, generally speaking, the misalignment becomes larger for 
larger anisotropics. This is more clearly quantified in Figs. 12 and 13, where the angles 
with respect to the z-direction of the tangent vector to the string and of the force are shown 
as a function of the angle of the direction of motion. From Fig. 12 we see that the tangent 
vector to the string systematically 'lags behind' the direction of motion as the latter varies 
from being aligned with the z-direction to being aligned with the x-direction. Only in these 
two limits does the string profile align itself entirely with the velocity. Moreover, the larger 
the anisotropy the more the strings 'wants' to stay aligned with the z-direction, changing 
direction quickly only as 99 approaches 7r/2. From Fig. 13 we see that the behaviour of the 
force is similar, except that for sufficiently large anisotropies its direction does not vary 
monotonically with the direction of the velocity. 

5 Discussion 

We have analyzed the drag force exerted on a quark moving through an anisotropic A/" = 4 
super Yang-Mills plasma described by the metric (2.2). In this case the anisotropy is 
induced by a position-dependent theta term in the gauge theory, or equivalently by a 
position-dependent axion on the gravity side. One may therefore wonder how sensitive the 
conclusions may be to the specific source of the anisotropy. In this respect it is useful to 
note that the gravity calculation involves only the coupling of the string to the background 
metric. This means that any anisotropy that gives rise to a qualitatively similar metric 
will yield qualitatively similar results for the drag force irrespectively of the form of the 
rest of the supergravity fields. 

An example of a rather robust conclusion is the ultra-relativistic behaviour of the drag 
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Figure 6. Drag force as a function of the velocity for a quark moving through an anisotropic 
plasma with a/T = 12.2 (left column) and a/T = 86 (right column) along four different directions 
lying at angles (curves from top to bottom) ip — 0, 7r/6, 7r/3, 7r/2 with respect to the longitudinal 
direction z. F is plotted in the appropriate units to facilitate comparison with the isotropic result 
for a plasma at the same temperature (top row) or at the same entropy density (bottom row). The 
isotropic result is given in eqs. (3.27) and (3.28). 



force. We have seen that the anisotropic solution (2.2) yields a drag force that becomes 
arbitrarily larger than the isotropic one for all ultra-relativistic quarks except for those 
whose velocity is perfectly aligned with the transverse xy-plane. This follows from the fact 
that the near-boundary fall-off of the metric (2.2) takes the schematic form 

5m- = ^ (^M- + ^'^-^ + + •••)• (5.1) 

As V grows closer and closer to 1 the string worldsheet develops a horizon closer and 
closer to the AdS boundary at n = 0. As a consequence the physics in this limit is solely 
controlled by the near-boundary behaviour of the metric. Eq. (3.18) can then be solved 
using the asymptotic form (5.1) of the metric functions. Generically the solution to leading 
order is determined by the C'(n^)-terms and yields oc 1 - v^. Substituting in (3.22) one 
gets F c<. 1/(1 — v^), or equivalently F = fip with a momentum-dependent drag coefficient 
fi (X p. For example, we show in Appendix A that the metric (2.2) yields a drag coefficient 



.^Aa^cosV ON 
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Figure 7. Drag force as a function of the direction of motion (/s, measured with respect to the 
longitudinal direction z, for a quark moving through an anisotropic plasma with a/T — 12.2 (left 
column) and a/T = 86 (right column) at three different velocities (curves from top to bottom) 
V — 0.9, 0.7, 0.5. F is plotted in the appropriate units to facilitate comparison with the isotropic 
result for a plasma at the same temperature (top row) or at the same entropy density (bottom 
row). The isotropic result is given in eqs. (3.27) and (3.28). 

at large p, where M is the quark mass. In contrast, in the isotropic case of [28, 29] 
the 0(ii^)-terms in the metric are absent. This means that the solution of eq. (3.18) is 
ti^ oc 1 — f ^ and hence that in this case the drag force in the ultra-relativistic limit has 
a softer divergence -Fi^o oc l/Vl — v"^- Rewriting this in terms of the momentum gives 
F = fip with a momentum-independent constant in this case. For certain choices of the 
parameters (for example for ip = 7r/2 in our case) the C'(u^)-terms in eq. (3.18) may vanish, 
in which case F oc 1/ Vl — v'^. 

The above discussion makes it clear that the linear behaviour of the drag coefficient 
in the ultra-relavistic limit, ^ oc p, depends solely on two features of the solution: The 

(2) 

presence of the gl^! term in the near-boundary expansion of the metric, and the fact that 
this term be anisotropic. The latter condition is necessary because otherwise the metric 
at order would be boost-invariant and there would be no solution for Uc at this order. 
Note that adding temperature to an otherwise boost-invariant theory will only affect gjfj 

(2) (2) 

and thus is not enough to make gliil non-boost invariant, since gin! is only a function of 
the external sources to which the theory is coupled. 
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Figure 8. Drag force as a function of the anisotropy for a quark moving along the longitudinal 
z-direction, i.e. at ip — (left column) or along the transverse x-direction, i.e. at ip = tt/2 (right 
column), at four different velocities (curves from top to bottom) v — 0.9, 0.7, 0.5, 0.25. F and a are 
plotted in the appropriate units to facilitate comparison with the isotropic result for a plasma at 
the same temperature (top row) or at the same entropy density (bottom row) . The isotropic result 
is given in eqs. (3.27) and (3.28). 

Interesting solutions with non-zero gf^J include mass deformations of the A/" = 4 super 
Yang-Mills theory in which some of the adjoint fermions acquire a mass, i.e. mass deforma- 
tions associated to a dimension-3 operator. The mass term is dual to a scalar field on the 
gravity side whose presence induces 0(u^)-terms in the metric. A simple example is the 
so-called GPPZ flow [31], which describes a deformation of = 4 that preserves A/" = 1 
supersymmetry.^ A generic further deformation of the GPPZ induced by an anisotropic 
external source will give rise to an ultra-relativistic drag coefficient of the form fi ck p. 
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Figure 9. Drag force as a function of the anisotropy for a quark moving at w = 0.5 (left column) or 
at u = 0.9 (right column) along four different angles (curves from top to bottom) = 0, 7r/6, 7r/3, 7r/2 
with respect to the longitudinal direction z. F and a are plotted in the appropriate units to facilitate 
comparison with the isotropic result for a plasma at the same temperature (top row) or at the same 
entropy density (bottom row). The isotropic result is given in eqs. (3.27) and (3.28). 
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A Ultra-relativistic limit 

In the limit v ^ 1 the value of Uc that solves eqn. (3.18) approaches the boundary, i.e. Uc 
0. Therefore in this limit Uc can be determined from of the near-boundary expansion of 
the metric functions, which takes the form: 

J- = 1 + ^ ( J-4 + ^ log n) + 0(n6) , 

B = l- + (^4 - ^ log.) + 0{u') , 

^ = 1 + - (— - ^ - ^ log.j + 0(n«) , (A.l) 



Figure 10. Projection of a t = snapshot of the string profile (3.7) (continuous curves) onto the 
a^z-plane for a quark moving with velocity u = 0.7 in four different directions (indicated by the 
dashed straight lines) that lie at angles (clockwise) f = 7r/18, 7r/6, 7r/3, Stt/IS with respect to the 
z-direction. The quark moves through a plasma with anisotropy a = 12.2 T (left) and a — 86T 
(right). The origin {x,z) = (0,0) corresponds to the string endpoint, which lies at the boundary 
M = 0. The coordinate u increases along the curves away from this point. 




0.05 0.10 0.15 0.20 0.25 0.30 ^ 0.05 0.10 0.15 0.20 0.25 0.30 ^ 

Figure 11. Generic misalignment between the direction of the quark velocity (dashed straight 
lines), the direction of the force (arrows) and the direction tangent to the string profile at its endpoint 
(continuous straight lines). The quark velocity is v — 0.7, its direction lies at angles (clockwise) 
(p = 7r/18, 7r/6, 7r/3, 87r/18 with respect to the 2;-direction, and the anisotropy is a = 12.2 T (left) 
and a = 86T (right). 



The coefficients i34,J^4 depend on a,T and are related to the energy and the pressures of 
the plasma eqs. (35) in [20]. They are not determined by the near-boundary analysis but 
must instead be read off from a full bulk solution. 

We must solve (3.18) for Uc to leading order in 1 — v^. For cos 99 7^ we only need to 
consider the terms of 0{u'^) in the metric functions, and the solution is 

Substituting into (3.22) we obtain the drag force 
-' \/A v"^ 

F = cos^ (f ^ (sin 93, cos 99) . [cos 99 7^0] (A. 3) 

Sir 1 — 
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Figure 12. Tendency of the string to align itself with the longitudinal direction z for four different 
anisotropics (from top to bottom) a/T — 12.2,20.3,42.6,744. The angle lPt- is the angle between 
the z-axis and the tangent vector to the string at its endpoint, defined as in eqn. (3.23). The angle 
If is the angle between the z-axis and the velocity. The magnitude of the velocity is v = 0.7 (left) 
and V = 0.9 (right). 
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Figure 13. Correlation between the direction of the force, (fY^ and the direction of the veloc- 
ity, (p (with both angles measured with respect to the longitudinal direction z), for four different 
anisotropics (from top to bottom) a/T = 12.2,20.3,42.6,744. The magnitude of the velocity is 
V = 0.7 (left) and v = 0.9 (right). 



The divergence when f — t- 1 contrasts with the softer behaviour (3.27)-(3.28) of the isotropic 
case. We conclude that for cosip ^ the ratio F^nisoZ-Piso diverges in the limit u — )■ 1 as 
l/Vl — w^, in agreement with our numerical results displayed in Fig. 6. Note that this is 
true even if the two plasmas have different temperatures and/or different entropy densities, 
since in the anisotropic case F diverges as 1/(1 — v'^) irrespectively of the temperature or 
the entropy density. 

The previous analysis shows that the limits u — ?■ 1 and (/J — >• 7r/2 do not commute. This 
is because if we first set cosy? = then the terms of order cancel out in eqn. (3.18) and 
we must go to order u'^. The solution in this case is 



,2 



1 r2 /121 Bi + J^i 

' [cos(/p = 0] (A.4) 



Vr^^ V 576 r4 T4 
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which yields the drag force 



2tT y/l 



121 
576 



rp4 



[cos (/J = 0] 



(A.5) 



This result is valid for any value of a/T, large or small, and it implies that the ratio 
-PlnisoZ-^iso is finite in the limit v — )• 1 and given by 



F 



1 



121a4 



7r2 \ 576r4 



T4 + B4 

rp4 



[cos (p = 0] 



(A.6) 



This result is valid for any a, large or small (as long as the motion is exactly aligned 
with the x-direction) . In order to proceed further we need analytic expressions for the 
coefficients J- 4,84. These are known in the limiting cases of small and large a/T. In the 
first case they are given in eqn. (175) of [20]: 



^4 

B4 



16 



16 



593 
1152 



"101 




-llo 


384 


12 


7 , 


/27rr^ 


— log 




12 ^ 


V a J 



27rr 

a 



7 /a 
— log — 
12 ^ VA 



+ 0(a6) , 
+ 0{a^) , (A.7) 



where A is a reference scale related to the conformal anomaly. Substituting into (A.6) we 
find 



Ik 



1 + 



+ 



[cos f = 0, small a/T] 



Note that the dependence on the reference scale A has cancelled out in this result, as 
expected from the discussion in the last paragraph of Sec. 2. The result (A. 8) shows 
that the drag force on an ultra-relativistic quark moving along the transverse directions 
in an anisotropic plasma with small a/T is greater than the drag in an isotropic plasma 
at the same temperature, in agreement with our numerical results. In order to make this 
comparison at equal entropy densities we use the fact that the entropy density at small 
a/T is given by (see eq. (174) in [20]) 



^ + " 



16 



a^ + 



Inverting this relation. 



1/3 



a/3 



1 

24 



NT] 
2tt J 



2/3 2 



+ 



substituting in (A.5) and taking the ratio with (3.^ 



s2/3 

we arrive at 



a 

s473 



(A.9) 



(A.IO) 



1 



1 

48 



iV2 

c 

2^ 



2/3 2 



32/3 



+ 



,4/3 



[COS If = 0,smana3/s] (A.ll) 



We see that, in contrast to the case of equal temperatures, the drag in the anisotropic 
plasma is smaller if the comparison is made at equal entropy densities, again in agreement 
with our numerical results. 
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In the limit of large a/T the coefficients J-4, can be obtained by combining eqs. (35), 
(89) and (90) of [20]. The result is 



-^4 



1 

132 
1 

6336 



132a^Ci„t + 77a'' log 



A 



-6336a^Ci„t + 1331a^ - 3696a* log 



348c,„t7r2a^/3r^^/3 + • 

4032c..,^2^i/3r"/3 



A 



+ 



(A.12) 
, (A.13) 



where Cint in an integration constant and Ccnt is the constant introduced in (2.5). Following 
the same procedure as in the small-a case we find that the ratio at equal temperatures is 



\/2^ 

rV3 



+ 



vr 



[cos 99 = 0, large a/T] 



(A.14) 



where the dots stand for subleading terms in the large a/T limit, and at equal entropy 
densities it is 



Ik 

F 



1 



1/48 



1 



16 



[cos 93 = 0, large a^/s] (A. 15) 



We conclude that at large anisotropics the ultra-relativistic drag in the anisotropic case is 
always greater than the isotropic drag. 



B Small-anisotropy limit 



For small values of a/T analytic expressions for the metric functions can be found [19, 20] 
by perturbing around the isotropic case. The result is 

,4 



F{u) 



1 



u 



+ a'F2{u) + 0{a'') 



B{u) = l + a'^B2{u) + 0{a^), 

,21 



\og'H{u) 



a^u^ 



^log 



1 + 



u 



0{a^) , 



where 



B2(u) 



24ul 



Su^iul - u^) - lOn^ log 2 + 3ui + lu" log 1 + 



24 



lOn^ 



+ log 1 + 



(B.l) 
(B.2) 
(B.3) 

(B.4) 
(B.5) 



Using these expressions in the general formulas of Sec. 3 we obtain the correction to the 
isotropic result for the drag force at leading order in a/T . The result for the drag force 
along the longitudinal direction z is 



= Fiso(T) 



/ ^2 X 1 - ^2 + vi-v2 + (1 + ^;2) log + 

^^yf^) 247r2(l-t;2) 



+ 



2^4 



B.6) 
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whereas for the transverse direction x it is 

+ yjl - i;2 + (4^2 - 5) log f 1 + VT 



= FUT) 



1+U 



1 



247r2(l 



+ 



a 

ji4 



(B.7) 

The 0{a'^ /T"^) correction in (B.6) is positive for v £ [0, 1], whereas that in (B.7) is negative 
for Q < V < Vc and positive for < fc < 1, where Vc — 0.9. This means that, for small 
enough an anisotropy, the drag force along the longitudinal direction in the anisotropic 
plasma is always larger than the drag force in an isotropic plasma at the same temperature 
(but different entropy density). In the case of motion in the transverse direction the 
anisotropic drag is smaller than the isotropic drag for low v and larger for high v. This is 
in agreement with the numerical results of Sec. 4. 

In order to compare with an isotropic plasma at the same entropy density (but different 
temperature) we use the relation found in [20] for the entropy density of the anisotropic 
plasma: 



1 



87r2T2 



+ 



a 

ji4 



(B.^ 



Inverting this relation and substituting in (B.6) and (B.7) we get: 



FUs) 



F. 



1 + 



1 + 



24 

„2 



iV2 

27rs J 



(1 - t;2) + (1 + ^2) log(l + Vl -t^2) 



24 V27rs J 



1 



(1 - t;2) + (4?;2 



+ 



a 



5) log(l + Vl - ^2) 



^4/3 
4 



1 



+ 



a 



54/3 



(B.9) 



In this case the leading correction is positive for all v in z-direction and negative for all v 
in the x-direction. Again, this is in agreement with the numerical results of Sec. 4. 
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